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ISO/SWS observations of SN 1987A: II. A refined upper limit 
on the mass of 44 Ti in the ejecta of SN 1987A.* 
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Abstract. ISO/SWS observations of SN 1987A on day 3 425 show no emission in [Fe I] 24.05 fj,m and 
[Fe II] 25.99 fim down to the limits ~ 0.39 Jy and ~ 0.64 Jy, respectively. Assuming a homogeneous distri- 
bution of 44 Ti inside 2 000 km s" 1 and negligible dust cooling, we have made time dependent theoretical models 
to estimate an upper limit on the mass of ejected 44 Ti. Assessing various uncertainties of the model, and checking 
the late optical emission it predicts, we obtain an upper limit of ~ 1.1 x 10 -4 Mq. This is lower than in our 
previous estimate using other ISO data, and we compare our new result with other models for the late emission, 
as well as with expected yields from explosion models. We also show that steady-state models for the optical 
emission are likely to overestimate the mass of ejected 44 Ti. The low limit we find for the mass of ejected 44 Ti 
could be higher if dust cooling is important. A direct check on this is provided by the gamma-ray emission at 
1.157 Mev as a result of the radioactive decay of 44 Ti. 
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1. Introduction 

The late (t > 150 days) emission from supernova (SN) 
t ' 1987A probes the content of radioactive isotopes in the 
, supernova ejecta. Due to the proximity of the supernova, 
$H ' it has been possible to follow its evolution to very late 
. 5^ i stages when 44 Ti takes over from 56 Co and 57 Co as the 
most important radioactive isotope for the powering of 
the light curve, i.e., after ~ 1 500 — 2 000 days (Woosley et 
al. 1989; Kumagai et al. 1991; Kozma & Fransson 1998a; 
Kozma 2000, henceforth K00). Both 56 Co and 57 Co are 
decay products of more short-lived nickel isotopes with 
the same mass numbers. Measuring the radioactive input 
therefore provides estimates of the ejected masses of 56 Ni, 
57 Ni and 44 Ti, i.e., Af( 56 Ni), M( 57 Ni) and M( 44 Ti), re- 
spectively, which constrain models of the explosion (e.g., 
Woosley & Weaver 1995; Thielemann et al. 1996; Nagataki 
et al. 1997). 

While M( 56 Ni) and M( 57 Ni) are rather accurately 
known for SN 1987A, with M( 56 Ni) w 0.071 M (e.g., 
Suntzeff & Bouchet 1990), and M( 57 Ni) - 3.3 x 10~ 3 M 
(Fransson & Kozma 1993, and references therein), the sit- 
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uation is more uncertain for 44 Ti. To some extent this has 
been due to a poorly known decay time of this isotope, 
although this has recently improved; 44 Ti decays first to 
44 Sc, on a time scale of 87.0 ±1.9 years (Ahmad et al. 
1998; Gorres et al. 1998; see also Norman et al. 1998 who 
obtain 86 ± 3 years), and then quickly further to 44 Ca. 

Based on time dependent modeling, K00 estimates 
that M( 44 Ti) = (1.5 ± 1.0) x 10~ 4 M© best explains 
the evolution of the optical broad-band photometry of SN 
1987A for t < 3 270 days. This is consistent with Chugai 
et al. (1997), who find Af( 44 Ti) - (1 - 2) x 10~ 4 M from 
the optical line emission at 2 875 days, and Mochizuki 
& Kumagai (1998) who obtain the same result from a 
light curve analysis (sec also Nagataki 2000). However, as 
is shown by K00, both broad-band photometry and the 
emission modeled by Chugai et al. (1997) constitute only 
a minor fraction of the total emission put out by the su- 
pernova. At these epochs most of the emission from the 
supernova instead comes out in the far infrared (IR) in a 
few iron lines, most notably [Fe II] 25.99 /im. This makes 
bolometric corrections to the late optical data rather un- 
certain, and a more direct way to measure the content of 
M( 44 Ti) is to measure the flux in the far-IR lines. 

In a recent study Lundqvist et al. (1999; henceforth 
L99) analyze Infrared Space Observatory (ISO; Kessler 
et al. 1996) data obtained mainly at t = 3 999 days. 
L99 couple these observations to time dependent model 
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calculations similar to those in K00, and from the ab- 
sence of iron line emission in the ISO spectra they de- 
rive M( 44 Ti) < 1.5 x 1(T 4 M . This limit is consistent 
with the results of Chugai et al. (1997), Mochizuki & 
Kumagai (1998) and K00. A different limit was obtained 
by Borkowski et al. (1997) who in their preliminary analy- 
sis of ISO data from t = 3 425 days derive M( 44 Ti) < 1.5 x 
10~ 5 M Q . Such a low titanium mass would have impor- 
tant consequences for models of the explosion. Here we 
analyze the data of Borkowski et al. in the same way as 
was done for the 3 999 day data in L99, however, expand- 
ing the discussion on the uncertainties of our method. In 
particular, we discuss the uncertainty of the temperature, 
which is important since the far-IR lines are collisionally 
excited, and their excitation energy (e.g., T « 550 K for 
the 26 /im line) is much higher than the gas temperature. 
Because the emission at this epoch is powered mainly by 
positrons from the decay of 44 Ti, we also discuss the un- 
certainty of the line fluxes due to the efficiency of the 
trapping of positrons. 

2. Observations and Results 

We have used the ISO data discussed in Borkowski et al. 
(1997). These data were obtained on 10 July, 1996 (i.e., on 
day 3 425 of SN 1987A) in the SWS06 mode of the Short 
Wavelength Spectrograph (SWS; de Graauw et al. 1996). 
The data we analyze span the regions 23.59 — 25.04 /im 
and 25.51 — 26.53 fj,m. Both these ranges are covered by 
band 3D of SWS, and include the lines [Fe I] 24.05 /im 
and [Fe II] 25.99 /im, i.e., the two strongest lines expected 
from the supernova (cf. L99). These data are superior to 
those discussed in L99 as the SWS06 scans are deeper 
than the SWS01 scans in L99, and include a large enough 
wavelength interval for accurate continuum determination 
compared to the SWS02 data in L99. 

The data were retrieved from the ISO archive, 
and have been reduced with the off-line processing 
(OLP/pipeline) version 8.7. We used ISO Spectral 
Analysis Package (ISAP) to flatfield and average over the 
scan directions. The averaging over the detectors was done 
with a resolution of R = 300, using the standard 2.5cr- 
clipping. 

In Fig. 1 we present fully reduced spectra of SN 1987A 
for the two wavelength regions. Although the nominal in- 
strumental resolution of the SWS spectra is R w 1 800 
we have averaged the spectra with a bin size of R = 
300, corresponding to w 1000 km s . This is sufficient 
to resolve the lines since they could extend to well over 
2 000 km s _1 . However, we detected neither [Fe I] 24.05 /im 
nor [Fe II] 25.99 /im in the spectra. We estimated the RMS 
in the data with a zero-order baseline fit in the regions 
23.80 - 24.80 /mi and 25.60 - 26.40 /im, and found 0.11 Jy 
and 0.19 Jy, respectively. 

The exact RMS-values vary somewhat with the bin size 
used in the averaging routine. Reducing the bin size to R 
= 500 gives the limits 0.13 Jy and 0.21 Jy, respectively. 
In order to derive a conservative limit on the 44 Ti mass 
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Fig. 1. ISO SWS/AOT6 spectra of SN 1987A on day 
3 425. The spectra were reduced using the pipeline soft- 
ware (see text). The bin size is w 0.081 (0.086) /im for 
[Fe I] 24.05 /im ([Fe II] 25.99 /im) (see text), which is low 
enough compared to the expected line width. Overlaid 
(dashed lines) are the line emissions from a model with 
M( 44 Ti) = 5 x 10~ 5 Mq (without photoionization) de- 
scribed in Table 1 and Sect. 3.2. To obtain the modeled 
line profiles we have assumed a maximum core velocity 
of 2 000 km s _1 , and a homogeneously distributed emis- 
sion throughout the core (see Sect. 3.2). The peak flux of 
[Fe II] 25.99 /im for this model is slightly less than what 
is needed for a 3a detection (cf. Sect. 3.3). 



we have adopted these values. A 3a limit for the 26 /im 
line then becomes ~ 0.64 Jy, while for the 24 /im line it is 
lower, ~ 0.39 Jy. The bin size is small enough so that we 
can use the limits for the peak of the expected line profiles 
in our estimate of the Ti mass. 



3. Interpreting the SWS observations 

3.1. The model for line emission from the supernova 

To interpret the upper limits on the iron line fluxes in 
Sect. 2, we have made model calculations similar to those 
in L99. Our computer code is described in detail in Kozma 
& Fransson (1998a), and summarized in L99. Here we only 
give a very brief recapitulation of the model, and discuss 
recent improvements. 
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The thermal and ionization balances are solved time- 
dependently, as are also the level populations of the most 
important ions. A total of ~ 6 400 lines arc included in 
the calculations. 

The radioactive isoptopes included are 56 Ni, 57 Ni, and 
44 Ti, and we calculate the energy deposition of gamma- 
rays and positrons solving the Spencer-Fano equation (see 
Kozma & Fransson 1992). We assume that the positrons 
deposit their energy locally. 

We include 0.07 M of 56 Ni and 3.33 x 10~ 3 M Q of 
57 Ni (Sect. 1). In L99 three values of M( 44 Ti) were tested: 
M( 44 Ti) = 0.5 x 10~ 4 , 1.0 x 10~ 4 and 2.0 x 10~ 4 M Q . Here 
we extend this to include three new values: M( 44 Ti) = 
1.0 x 10~ 5 , 1.7 x 10~ 5 and 3.0 x 10~ 5 M . Although wc 
follow the complete evolution of the supernova after 150 
days, we concentrate our discussion mainly on 3 425 and 
3 999 days, i.e., the epochs of the ISO observations. 

In L99 an e-folding time of 78 years was used for 
44 Ti, as this was thought to represent a mean value from 
experiments. More accurate measurements of the decay 
time (Ahmad et al. 1998; Gorres et al. 1998; Norman et 
al. 1998) became available while the analysis of L99 was 
completed, and the results of L99 were corrected accord- 
ingly before publishing. Here, we use the more accurate 
c-folding time (87.0 ±1.9 years) from the outset. 

The explosion model we use is the same as in L99 and 
Kozma & Fransson (1998a). That is, we take the abun- 
dances from the 10H model (Woosley & Weaver 1986; 
Woosley 1988), but distribute the elements so that hy- 
drogen is mixed into the core. Spherically symmetric ge- 
ometry is assumed, and the iron-rich core extends out to 
2 000 km s _1 , outside of which (out to 6 000 km s _1 ) a 
hydrogen envelope is attached. 

We use the Sobolev approximation for the line trans- 
fer. This is a good approximation for well-separated lines 
in an expanding medium, but is poorer when lines overlap. 
This is actually the case in the UV, especially at earlier 
epochs. The overlap leads to UV-scattering which affects 
the UV-field within the ejecta. As in L99, we have studied 
two extreme cases to test this effect: photoionization as in 
the original model, and simply switching off the photoion- 
ization caused by the UV-field. From the results of L99 we, 
however, do not expect the UV-field to be the dominant 
source of uncertainty in our models. Another drawback of 
the Sobolev approximation is that we do not account for 
line fluorescence in which UV lines are split into optical 
and IR lines. To treat this accurately is beyond the scope 
of this paper, but is treated in detail in a forthcoming 
paper. 

Since L99 we have updated the collision strengths for 
the [Fe II] 26 /im transition. The collision strength, O, is 
given by Zhang & Pradhan (1995) as a function of temper- 
ature. The value of Q at T < 200 K is set to ft = 7.0 (A. 
Pradhan, private communication). We note, however, that 
there can still be uncertainty in this result as a fully rcl- 
ativistic calculation is needed to accurately calculate the 
collision strength at these low temperatures (M. Bautista, 
private communication) . 
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"Distance = 50 kpc. Homogeneous distribution of emitting 
gas inside ±2 000 km s _1 has been assumed. The resulting 
shape of the line profile is seen Fig. 1. The flux values are the 
peak values of the lines assuming this profile, 
'indicates whether or not photoionization has been included 
in the calculation. (See Sect. 3.1.) 



3.2. Model calculations 

The results of our calculations can be seen in Table 1. We 
tabulate the fluxes at 3 425 days for [Fe I] 24 /xm and [Fe II] 
26 fim for the three values of M( 44 Ti) discussed in L99: 
5 x 10" 5 , 1 x 10~ 4 and 2 x 10~ 4 M©. We also include the 
three new values: 1 x 10" 5 , 1.7 x 10" 5 and 3 x 10" 5 M©. 
For each value of M( 44 Ti) we tabulate results for models 
with and without photoionization. In all models a simple 
form of dust absorption was assumed (Kozma & Fransson 
1998b). 

The line fluxes in Table 1 are for a distance to the su- 
pernova of 50 kpc, and assuming a homogeneous distribu- 
tion of emitting gas inside ±2 000 km s _1 . The FWHM 
velocity for the resulting profile (shown in Fig. 1) is 
»fwhm = 2 828 km s _1 , which is close to what Haas et al. 
(1990) observed for [Fe II] 17.94 fim at - 400 days after 
the explosion, ufwhm = (2 900 ± 80) km s _1 . 

In L99 it was discussed how the fluxes in [Fe I] 
24.05 (im and [Fe II] 25.99 /xm (/W^m and /26/xm, re- 
spectively) scale with increasing values for M( 44 Ti) at a 
given epoch. It was found that /26^m scales nearly linearly 
with M( 44 Ti), while /24 M m has a weaker dependence. The 
stronger dependence for /26/jm is due to an ionization effect 
because A^eii, the relative fraction of iron in Fe II in the 
iron-rich gas, increases with increasing M( 44 Ti). We find 
that this trend continues also for M( 44 Ti) < 5 x 10~ 5 M©. 
This is shown in Fig. 2 and 3, where models joined by 
dashed (solid) lines are without (with) photoionization. 
Whereas Xpeii and the temperature in the iron-rich gas 
are only shown for 3425 days (Fig. 3), results for the line 
fluxes (Fig. 2) are shown for both 3425 and 3 999 days. 
Note how the fluxes remain nearly constant between these 
two epochs, especially in the case of /26/xm- 
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Fig. 2. Modeled fluxes of [Fe I] 24.05 and [Fe II] 25.99 fim from SN 1987A as functions of the mass of ejected 44 Ti, 
3 425 and 3 999 days after explosion. Models with photoionization included (see text) are shown by filled squares and 
joined by solid lines, while for the others photoionization was not included. The inferred 3<r detection limits for the 
two epochs are shown as dashed-dotted lines. Note the near constancy of the modeled fluxes between the two epochs, 
and the rather small uncertainty due to photoionization. 



The only deviation from a linear log(M( 44 Ti)) ver- 
sus log(/) behavior in Fig. 2 is for /24 Mm at the high- 
est values of M( 44 Ti) in the models with photoioniza- 
tion. This is because here Xp c i (the fraction of iron in 
Fe I in the iron-rich gas) starts to fall significantly below 
unity (cf. Fig. 3): for M( 44 Ti) = 0.3, (0.5), (1.0), (2.0) x 
10~ 4 M , X Fcl « 0.78, (0.70), (0.58), (0.44) at 3 425 days. 
(The same number for the case without photoionization 
is X Fd « 0.82, [0.76], [0.68], [0.58].) The temperature in- 
creases monotonically with M( 44 Ti) (Fig. 3) and affects 
/24/xm and /26/jm in the same way. This is because the two 
lines have nearly the same excitation energy, and their 
effective collision strengths vary only weakly (and in our 
models are assumed to be constant over the small temper- 



ature regime [~ 115 — 180 K] found in the iron-rich gas in 
our models). 

Table 1 and Fig. 2 show that the 26 /an line is stronger 
than the 24 jum line even for the models with the lowest 
values of M( 44 Ti) (and thus the highest values of Xpd). 
This is because the collision strength for [Fe II] 25.99 //m 
line is much larger than for [Fe I] 24.05 /im. 

Table 1 and Fig. 2 also show that switching off pho- 
toionization does not have a dramatic effect on the line 
fluxes. Again, the largest difference for the fluxes is for 
the 24 fim line and at the highest M( 44 Ti) considered. 
For values < 5.0 x 10 -5 M©, switching off photoionization 
only affects / 24([1 m (J^m) by < 11% ( < 22%) at 3 425 
days. The small difference is simply due to the minor shift 
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Fig. 3. Fraction of iron in Fe II, Xpeii, and temperature in the iron-rich part of the ejecta of SN 1987A at 3 425 days 
as a function of mass of Ti. Symbols and lines have the same meaning as in Fig. 2. As expected, the temperature 
and ATpeii increase monotonically with increasing mass of 44 Ti. 



in the degree of ionization, from Fe II to Fe I, when pho- 
toionization is switched off. 

3.3. Mass of 44 Ti from [Fe II] 25.99 (im 

With the results in Fig. 2 it is straightforward to estimate 
the upper limit on M( 44 Ti). It is clear that the best es- 
timate within the framework of our modeling comes from 
the 26 /im line on day 3 425. Using ~ 0.64 Jy from 

Sect. 2.1, we find the upper limit M( 44 Ti) < 5.2 (5.9) x 
10~ 5 M Q with (without) photoionization included. A con- 
servative limit (i.e., the case when photoionization is 
unimportant) is therefore M( 44 Ti) < 5.9 x 1CT 5 M©. In 
Fig. 1 we have included the expected line emission for a 
model with M( 44 Ti) = 5.0 x 1(T 5 M , i.e., close to this 
limiting mass. Although our limit is nearly a factor of four 
higher than that found by Borkowski et al. (1997), using 
the same data, it is still substantially lower than found by 
L99 from day 3 999. Using our updated code, the day 3 999 
data only set a limit M( 44 Ti) < 1.1 x 1CT 4 M©. This is a 
factor ~ 1.7 better than can be obtained from the 24 /im 
line on day 3 425. We will evaluate our limit from the 
26 /zm line on day 3 425 in Sect. 4.1. 

4. Discussion 

4.1. Uncertainties in the modeling 

There are several uncertainties involved in our modeling 
of the line fluxes. The effect of photoionization was al- 
ready discussed in the previous section, and was found 
to be rather small, as was also argued for in L99. An 
even smaller uncertainty in the estimated M( 44 Ti) is in- 
troduced by the decay time of 44 Ti. Given the small error 
in this (Ahmad et al. 1998; Gorres et al. 1998; Norman et 



al. 1998), we can safely ignore this uncertainty. More un- 
certain is the distance to the supernova. This is uncertain 
to - 10% (Walker 1998; Lundqvist & Sonneborn 2001), 
and a conservative limit is given by Gibson (2000), cor- 
responding to an uncertainty in the line flux of ~ 27%. 
Because the modeled ./26/jm scales almost linearly with 
M( 44 Ti) (Sect. 3.2), we have used this number also for 
the uncertainty in the derived M( 44 Ti). 

A source of uncertainty, not investigated in L99, is 
the rather unknown rates of charge transfer for the ele- 
ments in the core. In the models in Sect. 3, no charge 
transfer was included between Fe and excited states of 
He, as well as between atoms and ions of Fe. We have 
tested the importance of including this charge transfer, 
using the rates suggested by Liu et al. (1998) in a model 
with M( 44 Ti) - 1.0 x 10~ 4 M and no photoionization. 
Compared to the model with the same M( 44 Ti) in Sect. 
3, the model with the full charge transfer has a higher 
value of XfcII by ~ 20%, mainly due to charge transfer 
between Fe I and Fe III. However, the temperature in the 
model with full charge transfer is somewhat lower, so the 
difference in /26/im between the two models is only ~ 7%. 
With full charge transfer included, the estimated M( 44 Ti) 
would therefore be lower than in Sect. 3, but due to the 
uncertainty of the charge transfer rates, and in order to 
derive a conservative limit on M( 44 Ti), we assign a gen- 
erous error of 10% in M( 44 Ti) due to this effect. 

L99 studied the uncertainty of the collision strength of 
the 26 fim line. Although we are now using a better fit 
to the results of Zhang & Pradhan (1995), the collision 
strength at very low temperatures is still uncertain (see 
Sect. 3.1). In L99 we assumed an uncertainty of 15% in 
M( 44 Ti) due to this, and we retain this number also here. 

The collisional excitation of the 24 /im and 26 /im 
lines is more sensitive to the temperature (through the 
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exp(— hc/XkT) term) than to the collision strength. Too 
high a temperature in our models could therefore over- 
estimate the flux in these lines, and we would conse- 
quently underestimate M( 44 Ti). We can check this ef- 
fect by studying Table 1 where we see that models with 
M( 44 Ti) = 5.0 x 1CT 5 M and M( 44 Ti) = 1.0 x 10~ 4 M Q 
differ in /26^m by a factor of ~ 2.2 at 3 425 days. Because 
the temperature in the M( 44 Ti) = 1.0 x 10~ 4 M model 
is w 159 K (see Fig. 3), a lowering of the temperature by 
only ~ 20% would decrease to the same level as in 

the M( 44 Ti) = 5.0 x 10" 5 M Q model, for constant X FcU . 
However, the temperatures in our models are not free pa- 
rameters but are fixed by the heating and cooling. The 
energy not coming out in the 24 /an and 26 //m lines must 
instead come out in other lines (like the optical/IR recom- 
bination lines) which are less sensitive to temperature, or 
as emission from dust. We will make a consistency check 
of the optical and far-IR emission in our models in Sect. 
4.2, but we note already here that only a small transfer 
of energy loss from the far-IR lines to the optical/IR re- 
combination lines would increase the optical and IR flux 
considerably, and then the models of K00 for the optical 
would severely overestimate M( 44 Ti). We therefore believe 
we are not making a serious error in the derivation of the 
temperature in our models. 

Dust emits a continuum which is difficult to detect, 
but we may for simplicity assume that the dust emission 
does not interact with the gas. In that case, dust cooling 
has exactly the same effect for the 24 /im and 26 /im line 
fluxes as just lowering M( 44 Ti). (See the models in Sect. 
3). It may therefore well be that dust cooling could tap 
the supernova of its [Fe II] 25.99 ^m emission, but it is 
likely that it would do so in a way which would also make 
the optical and near-IR lines too weak. We do not assign 
an explicit error in the estimated M( 44 Ti) due to dust 
cooling. (See also Sect. 4.2.) 

L99 discussed uncertainties in the modeled line fluxes 
due to the explosion model used. As in L99 we have 
used the 10H explosion model (Woosley & Weaver 1986; 
Woosley 1988), mixed by Kozma & Fransson (1998ab) to 
give good agreement between their modeling and late time 
observations. Kozma & Fransson (1998b) compared the 
results from this model with a similarly mixed version of 
the 11E1 model (Shigeyama et al. 1988), and found that 
the iron lines are insensitive to the explosion model used. 
This is because the iron core mass, which is fixed by the 
amount of ejected 56 Ni, is the same in both models. We 
follow L99 and assign a 15% error in M( 44 Ti) due to the 
choice of explosion model. 

In our calculations we have assumed a local deposition 
of positron energy from the radioactive decay of 44 Ti. This 
assumption was also used in L99. The argument is that 
the efficiency of trapping cannot change substantially until 
day ~ 3 600 (K00), and that the most straightforward in- 
terpretation of a near-constant efficiency of trapping over 
an extended period of time is to assume full trapping. L99 
assumed a 15% error in M( 44 Ti) due to the uncertainty of 
trapping. Here we do not infer an explicit error to the mod- 



els in Sect. 3, but in Sect. 4.2 we investigate this in greater 
detail in terms of a consistency check of the modeled IR 
and optical fluxes. The positron deposition is discussed 
in Sect. 4.2 in conjunction with clumping. Even without 
positron leakage, clumping of the iron-rich gas could cause 
an error in the estimate of M( 44 Ti). However, according 
to L99, this error is likely to be small, ~ 5%, and can be 
ignored when compared to other errors discussed above. 

The combined error of f26fim due to model approxima- 
tions, except for different clumping and positron leakage 
scenarios (which will be discussed in Sect. 4.2) is therefore 
~ 36%. With the line profile used in Fig. 2 we thus arrive 
at an upper limit on M( 44 Ti) which is - 8 x 10~ 5 M Q . 

4.2. Consistency check of the modeled optical and 
infrared emission 

Although the limit on M( 44 Ti) we found above in Sect. 
4.1 is a factor of ~ 2 lower than the limit found by L99 for 
the 3 999 data, it is still compatible with K00 (see Sect. 1). 
However, the agreement is not perfect although the same 
computer code has been used. Here we make a consistency 
check of the modeled optical and IR emission produced by 
our model. For the late optical photometry (day 3 268) we 
adopt the same data as used by K00 (Soderberg, Challis & 
Suntzeff 1999). The observational error of these is approx- 
imately ±0.05 magnitudes, as judged from the system- 
atic error in the zero point of HST/WFPC2 photometry 
(Mould et al. 2000; P. Challis 2000, private communca- 
tion). 

We have concentrated our consistency check on models 
with M( 44 Ti) = 10~ 4 M . This is close to the upper limit 
found in Sect. 4.1 from the far-IR lines when the errors 
have been considered. This value of M( 44 Ti) is also in 
the lower range of the preferred values found by K00 in 
her modeling of the optical and near-IR emission. This 
can be seen from Table 2 where the modeled V, R and / 
magnitudes in model Ml arc all slightly fainter, by roughly 
half a magnitude, than the observed, while the predicted 
[Fe II] 26 /im line emission is somewhat too strong to be 
accommodated by our error analysis in Sect. 4.1. (Model 
Ml is identical to the model with M( 44 Ti) = 10~ 4 M Q 
and no photoionization in Table 1.) 

It is obvious that our models contain simplifications 
that can make us systematically overproduce far-IR emis- 
sion at the expense of optical emission. We still have to ex- 
plore fluorescence (see Sect. 3.1), and other simplifications 
could be clumping and different scenarios for the deposi- 
tion of positron energy. These two effects are probably also 
coupled. In Sect. 4.1 we mentioned that L99 found that the 
estimated error in M( 44 Ti) due to clumping alone should 
be small. This is true if one looks at the total emission 
coming out from the core, and if the energy deposition 
is local. However, if we assume that the positron optical 
depth scales linearly with density (ignoring effects due to 
the magnetic field) , positron leakage from the low-density 
component is more likely than from the high-density one. 



P. Lundqvist et al.: Mass of 44 Ti in the ejecta of SN 1987A 



7 



To test this we have, similar to what we did in L99, used 
a two-component model with equal mass in the dense and 
the less dense components. For the density ratio we use the 
value 5. From a purely geometrical point of view, positrons 
produced in a high-density clump see an optical depth for 
absorption in the positron emitting clump which is ~ 3 
times higher than for self-absorption of positrons emit- 
ted in a low-density region. It is therefore more likely that 
positrons from the low-density gas leak from their produc- 
tion site, and the chance that they are captured in either 
a low-density, or high-density component on a larger scale 
is about equal. From a physical point of view, absorption 
of the leaking positrons is slightly biased toward absorp- 
tion in the low-density gas since this is more ionized (due 
to less efficient recombination), which results in slightly 
higher deposition there. We have tested three scenarios, all 
without photoionization and with M( 44 Ti) = 1CP 4 M : 
the first (M2 and M3) is the same as in L99 with local de- 
position of the positron energy. The second (M4 and M5) 
is with all positron deposition (i.e., from both the low 
and high-density gas) occurring only in the high-density 
clumps. This situation is an extreme case (cf. above), but 
could be possible if magnetic effects make deposition in 
the high-density gas more likely than in the low-density 
gas. The third (M6 and M7) scenario is for a situation 
where the positrons in the low-density gas are instead as- 
sumed to leak into and deposit their energy in the the 
silicon-rich gas which is macroscopically mixed with the 
iron-rich gas. We have run these models both as steady 
state and time dependently to separate this effect as well. 
We have summarized the results in Table 2, where we show 
/24/xm, /26/^m, V, R and I at 3 425 days. (The observed op- 
tical fluxes in Table 2 are for 3 268 days, and not for 3 425 
days. The change in modeled optical flux between these 
two epochs is, however, negligible.) 

From the results in Table 2 it is clear that the far- 
IR emission is rather robust to density inhomogeneities 
in the core, as well as to the assumption of steady state 
as long as the positrons are all trapped in the iron-rich 
gas. Allowing for positron diffusion into the Si-rich gas 
decreases the [Fe II] 25.99 /xm emission; instead the emis- 
sion in V increases notably, while R and / are less af- 
fected. In our most extreme model M6 (cf. Table 2) the 
positron escape leads to a reduction of /26^m by ~ 32%. 
This model actually produces 26 fim emission which is be- 
low the detection threshold, so this shows that for extreme 
situations of clumping and positron depositions our model 
could accommodate up to M( 44 Ti) ~ 1.1 x 10~ 4 M Q in 
SN 1987A, and still not produce detectable far-IR emis- 
sion. However, such a scenario would imply V ~ 19.5 at 
3 425 days, a factor ~ 1.6 higher flux in V than observed. 
As the optical flux is only a minor contribution to the to- 
tal emission coming out from the supernova at this epoch, 
and we yet have to investigate fluorescence to model the 
optical emission in detail, we cannot exclude a scenario 
like this. We therefore base our upper limit on the 26 /im 
line which we believe we have modeled more accurately. 
To estimate a conservative upper limit on M( 44 Ti), taking 



Table 2. Modeled luminosities and magnitudes at 
3 425 days a 



Model Time dep. 6 / 24flm c f 2G „m c V R I 



Ml d 


yes 


0.32 


1.18 


20.4 


19.7 


20.0 


M2 e 


yes 


0.39 


1.30 


20.3 


19.7 


20.0 


M3 e 


no 


0.35 


1.27 


20.3 


20.4 


20.2 


M4 / 


yes 


0.47 


1.13 


20.1 


19.6 


19.9 


M5 / 


no 


0.46 


1.15 


20.1 


20.4 


20.2 


M6 9 


yes 


0.32 


0.80 


19.6 


19.5 


19.7 


M7 S 


no 


0.30 


0.82 


19.6 


20.0 


19.9 



a M( 44 Ti) = 10~ 4 M Q . No photoionization. Observed magni- 
tudes (at 3 268 days) are V = 20.0, R = 18.9, and I = 19.3 
(Soderberg et al. 1999, P. Challis 2000 [private communca- 
tion]). 

b Time dependent or steady state model. 

c Peak flux in Jy for the line profile used in Figs. 1 and 2. 

Observed upper limits on /24pm and /24/jm are 0.39 Jy and 

0.64 Jy, respectively (see Sect. 1). 

d One-component model, i.e., same as in Table 1. 

e Two-component model. Positron deposition in both density 

components. 

^Two-component model. Positron deposition only in the 
high-density component. 

9 Two-component model. Positron deposition in the high- 
density component and in the Si-rich gas (see text). 



into account various clumping and deposition effects, we 
thus arrive at an upper limit of 1.1 x 10~ 4 M©. 

4.3. Considerations of the derived mass of 44: Ti 

Models for the yield of 44 Ti give quite different results. 
In particular, in the model of Timmes et al. (1996; see 
also Woosley & Weaver 1995) with a zero-age mass of 
Mzams = 20 M (i.e., corresponding to SN 1987A) the 
mass of the initially ejected 44 Ti is 1.2x 10 -4 M©, but only 
1.4 x 1CT 5 M Q escape after fallback. This is in accord with 
our upper limit in Sect. 4.1, but it should be emphasized 
that the variation of M( 44 Ti) with Mzams in Timmes 
et al. (1996) is complex, and for models with Mzams = 
18 M and 22 M Q , the calculated M( 44 Ti) is closer to 
our limit. Furthermore, the recent models of Hoffman et 
al. (1999) for 15 and 25 solar mass models, indicate that 
the yield of 44 Ti in Woosley & Weaver (1995) may have 
to be increased. 

The models of Thielemann et al. (1996) have larger 
entropy and thus more alpha-rich freeze-out than those 
of Woosley & Weaver (1995). Accordingly, the ratio 
M( 44 Ti)/M( 56 Ni) (where M( 56 Ni) is the mass of ejected 
56 Ni that does not fall back) is higher. For example, in 
their 20 M Q model M( 56 Ni) w 0.074 M Q and M( 44 Ti) « 
1.7 x 10~ 4 Mq. This value for M( 44 Ti) is higher than our 
upper limit in Sect. 4.2. With the refinements made by 
Hoffman et al. (1999) to the models of Woosley & Weaver 
(1995), it seems that our upper limit is lower than, or at 
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least close to, M( 44 Ti) in the models of both groups. The 
calculated M( 44 Ti) is, however, very model dependent, 
and sensitive to, e.g., the explosion energy, which in the 
case of SN 1987A is still only known to an accuracy of 
- 30% (Blinnikov et al. 2000). 

The models of Nagataki et al. (1997, 1998) are related 
to those of Thielemann et al. (1996), although they also 
allow for 2-D. With no asymmetry, Nagataki et al. obtain 
M( 44 Ti) - 6 x 10~ 5 M for a SN 1987A-like explosion, 
which is a factor of ~ 3 lower than Thielemann et al. 
(1996), and could indicate the range of uncertainty of the 
modeling. In the models of Nagataki et al. the yield of 
44 Ti quickly increases with the degree of asymmetry. For 
an asymmetry of 2 between the equator and the poles, 
the stronger alpha-rich freeze-out in the polar direction 
increases M( 44 Ti) to M( 44 Ti) - (1.2-1.9) x 10~ 4 M Q (for 
M( 56 Ni) « 0.07 M Q ), the range of titanium yield depend- 
ing on the form of the mass cut (Nagataki 2000). Nagataki 
(2000) argues that this degree of asymmetry gives a dis- 
tribution of ejected 56 Ni which agrees with observed line 
profiles of [Fe II] 1.26 /xm (Spyromilio et al. 1990) and 
[Fc II] 18 /jm (Haas et al. 1990) at ~ 400 days. 

Although the models of Nagataki (2000) thus hint that 
asymmetry can explain the distribution of the ejected 56 Ni 
in SN 1987A, the need for asymmetry to reach 44 Ti in 
excess of 10~ 4 M Q cannot be considered a strong case, 
given the uncertainty in the modeling. We therefore do not 
regard our upper limit on M( 44 Ti) ~ 1.1 x 10~ 4 M Q for 
SN 1987A as contradictory when compared to the absolute 
yield of 44 Ti in, e.g., Nagataki (2000), especially since the 
models of Nagataki are trimmed to produce M( 44 Ti) in 
agreement with the light curve results of Mochizuki & 
Kumagai (1998). 

Table 2 also shows that although steady-state models 
can be used for the far-IR emission from the core, they fail 
to reproduce the optical broad-band emission. In particu- 
lar, they underproduce R and / magnitudes. For the cases 
in Table 2, the maximum errors are 0.8 and 0.3 magni- 
tudes, respectively. This is due to the freeze-out effect de- 
scribed by Fransson & Kozma (1993). Consequently, any 
estimate of M( 44 Ti) based on steady-state calculations for 
the optical broad bands is likely to overestimate M( 44 Ti), 
even though the structure and atomic data are correct. It 
could be that the models of Mochizuki & Kumagai (1998; 
see also Nagataki 2000) suffer from this, which may ex- 
plain why their lower limit, M( 44 Ti) = 1.0 x 10~ 4 M Q , al- 
most coincides with our upper limit, and why their lower 
limit is higher than the one of K00. 

Finally, we note that our upper limit on M( 44 Ti) could 
be higher if dust contributes significantly to the cooling at 
late times. This is not unique to our models, but affects all 
models trying to reproduce the optical and IR emission. 
The range of M( 44 Ti) bracketed by K00 and our analysis, 
M( 44 Ti) = (0.5 - 1.1) x 10~ 4 M would then be shifted 
to a range with higher masses. Heavy dust formation in 
the core would, however, block out the line emission from 
the core and affect line profiles. To test the importance 
of dust cooling, direct measurements of the gamma-ray 



emission at 1.157 MeV, the result of the radioactive de- 
cay of Ti, are needed. Future, and even more sensitive 
gamma-ray instruments than INTEGRAL, are needed for 
such a study. 

5. Conclusions 

We have extracted and modeled ISO/SWS data for day 
3 425 of SN 1987A. The supernova was neither detected in 
[Fe I] 24.05 /jm nor in [Fe II] 25.99 /im, which are the lines 
expected to be the strongest in the spectrum. For a likely 
line profile, we have used this finding to put an upper limit 
on the mass of ejected 44 Ti, M( 44 Ti). Including possi- 
ble modeling errors, we obtain M( 44 Ti) < 1.1 x 10~ 4 M Q . 
This supersedes the result of Lundqvist et al. (1999), who 
found a somewhat higher limit analyzing other data. The 
new limit is, however, more reliable in that we have in- 
cluded more possible sources of error, as well as late optical 
broad band data. It is shown that time dependent effects 
are important for the latter, probably explaining why our 
models give a somewhat lower preferred value of M( 44 Ti) 
compared to those in previous steady-state analyses. 
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